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Abstract. 

A parametric approach is developed to the method of S-tree diagrams 
and its generahzation for investigation of the hierarchical substructure of 
A^-body nonlinearly interacting systems (e.g., clusters of galaxies). The 
introduction of a parametric function allows us to take into account the 
individual features of the particles which do not have a direct influence 
on the dynamics of the system but are crucial for analyzing the output 
data. An algorithm is proposed for the latter based on the construction of 
parametrical subgroups. 

A generalization of the S-tree scheme is presented as well. This enables 
us to perform parallel substructure analysis with respect to various depen- 
dent or independent parameters, thus using fully the initial information of 
a system. 

1 Introduction 

Computers are powerful tools for the investigation and numerical simulation 
of complex systems, including information about their hierarchical struc- 
tures A class of such problems, namely N-body gravitating configura- 
tions, has an essential role in the understanding of various structures in the 
universe such as star clusters, galaxies, clusters of galaxies, and so forth. 

In the present paper we suggest an algorithm for a statistical method 
of investigating the hierarchical structure of gravitating N-body systems. 
The large-scale distribution of galaxies in the universe and the hierarchical 
structure of clusters of galaxies are two of the key problems of extragalactic 
astrophysics ||2|, ^. This is important because the relaxation times of clusters 
of galaxies are of the order or exceed the age of the universe, those systems 
recall their initial conditions. Hence, the study of the hierarchical properties 
of systems may directly indicate the mechanisms of their formation and 
evolution. 

Among the statistical methods used for investigating the hierarchical 
properties of galaxy clusters are the two-point correlation functions, Lee 
statistics, topological measures, wavelets and so on, for review see |0, |5[. 



The studies already performed enable us to conclude that many galaxy 
clusters may include subgroups with vivid dynamical peculiarities. The 
existence of the subgroups is also apparent from X-ray radiation data of the 
clusters. The development of more informative methods for the analysis of 
the problem is of particular interest, given the increase of the possibilities 
of observational technique. 

The S-tvee diagram method ^ is based on the idea of introducing 
a degree of boundness between the particles and their sets and therefore 
self-consistently using both positional and kinematic information about the 
particles. The idea is therefore that the parameters of particles in a phys- 
ically interacting system, such as the components of their 3-velocities and 
coordinates, have to correlate with each other. On the contrary, particles 
that do not satisfy the correlation have to be chance projections either in 
front of or behind the physical system. The problem is therefore to reveal the 
corresponding function of interaction, and to define its "degree" , since obvi- 
ously, various particles can have different degrees of interaction. This idea 
was realized via the S-tree technique, thus revealing the physically interact- 
ing system among the projected set of galaxies, and moreover establishing 
the subgroups with various degrees of interaction. This is achieved by con- 
sidering the geometrical properties of the phase space of the system, namely, 
its two-dimensional curvature. The N-body system as a hamiltonian system 
is thus being transferred via the Maupertuis principle to the phase space and 
its dynamics are being described by the flow of geodesies on a riemannian 
manifold with its metric determined by the potential of interaction . The 
S-tree method has already been used for the study of ESO Key Program 
data on the nearby Abell clusters and many other observational data. In 
particular, the ESO data enabled the determination of a new class of galaxy 
configurations as dynamical entities given their peculiar properties. 

The S-tree technique has been developed not only for point approxima- 
tion, but also to account for the extension of particles Q, and to include 
certain generalized schemes poj] . 

The parametric approach of the S-tree method developed in the present 
paper allows us to consider individual properties of the particles, which 
are crucial for the analysis of the output information on the hierarchical 
structure of the system. For example, in the case of the galactic systems 
such a quantity can be the morphological class, which though does not 
influence the dynamics directly, has an essential impact on analyzing the 
output. The aim is therefore to have the algorithm include the parametric 
representation in the hierarchic structure intrinsically. 
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Then we extend the approach to generahzed schemes that enable the 
algorithm to use different functions of boundness self-consistently, thus ob- 
taining a more complete picture about the hierarchical structure of the sys- 
tem. We first introduce the necessary definitions and then move on to the 
algorithm. 

2 S-tree and its generalization 

We now briefly recall the main steps of the S-tree method. As mentioned 
previously, the key idea of this method is the introduction of the concept of 
the degree of boundness p. 
Consider a set of N points: 

X = {xi, . . .,xn}, 

the function P 

P : X xX ^ R+ and p G R+. 

The basic definitions follow. 

Definition 2.1. We say that Vx G X and My ^ X are p-bounded, if 
P{x,y) > p. 

Definition 2.2. We say that U C X{U 7^ 0) is a p-bounded subgroup, if: 

1. Vx eU and Vy e U ^ P{x, y) < p; 

2. Vx G [/ and \/y & U 3x = Xi^,Xi^, . . . , Xi^ = y, 
that Xjj^j) > p; V/ = 1, . . . , A; — 1. 

Definition 2.3. We say that Ui, . . . ,Ud is the distribution of the set X via 
p-bounded groups, if: 

1. [jf=lU^=X■, 

2. i^j {i,j = l,...,d)^U,nUj=^; 

3. Ud{i = 1, . . . , d) is a p-bounded group. 

It is possible for the function P to consider different physical quantities. 
In various quantities such as the energy, the potential, perturbations of 
potential, momentum, and so on, are listed and discussed. For example, 
the choice of the mutual distance of the particles to define the subgrouping, 
obviously, is equivalent to the corresponding correlation functions (spatial 
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or angular, related with each other simply via the Limber equation), which 
is a rather incomplete characteristic for that aim. It can be shown that, 
at least for astrophysical problems, among the most informative ones is 
the riemannian curvature of the configuration space [Q], which determines 
the behavior of close geodesies, as known from basic courses on classical 



mechanics [11 



So, by the S-tree algorithm we obtained the distribution of the N-body 
system for any given function P and \/p. This distribution will satisfy Defi- 
nitions 2.1, 2.2 and 2.3. The final result can be represented either through 
tables or by graphs (S-tree). 

Now we aim to use the S-tree scheme for different boundness criteria 
simultaneously. 

We have t different functions Pi, . . . , Pt, which satisfy the definitions 
introduced in the method of S-diagrams and that describe the degree of 
interaction of the system's members. To each P^ we correspond a matrix 
where = {dfj),a = l,...,t^. 

Consider the following matrix D: 



D = {d^j), i,j = 1,...,A^, 

where 

dij = {d}j,dfj, d\j), i,j = 1,...,N, 

which contains the entire information of the system, that is, it uses all Pt 
functions. For example, the previously mentioned Riemannian curvature is 
a function of the coordinates, velocities of all the particles, their masses, 
potential of interaction, its derivatives, and so on. 

Each function Pq, a = 1, . . . ,t corresponds to a definite vector of bound- 
ness Pa = {pi, • • • , Pq^) and Qa satisfies the conditions of the choice of the 
vector of boundness ||5|. 

We construct the vector pD = {pi , ■ ■ ■ jPq^^), where = {pl^,. ■ ■ ,pij] 
k = l,...,QD; ka = l,.--,Qa, a = l,...,t. 

The vector p^i is obtained with the help of all possible variations of 
the components of p^{a = 1, . . . ,t;ka = 1, • • • , Qa), and it is obvious that 
Qd < {{N^ -N)/2 + lY. 

Taking into account the fact that different P^ functions contain different 
information about the system, we introduce the concept of the "degree of 
influence" Ea for each of the functions under consideration. We will con- 
ditionally understand as the number of physical quantities (distance, 
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potential of the interaction and its derivatives, and so on), existing in the 
Pa formulae. 

The next step of the algorithm is the transition from matrix D to matrix 
Du with the help of a fixed current vector p^. 



Du = {uij); i,j = l,...,N, 

"ij — V "-ij ) • • • ! "-i j ) i ''^tj ( 



where Uij = (uj^, . . . , u* •); uf Ja = 1, ... ,t) is defined in the following way: 



Ea if dS>pL 

if dg<pt- 

Thus, the matrix contains information on the degree of influence of 
each criterion too. Then we construct the matrix in the following way: 

= (vij); i,j = l,...,N, 

where Vij = YL=i ufj- 

The matrix characterizes the quantitative degree of boundness. At 
the same time we also take into account information on the p-boundness of 
particles, because matrix is defined with the help of matrix Z)„. Accord- 
ing to the S-tree method we then construct the vector of boundness fl for 
the matrix D^: 

p.= (^i,...,^t)- 

It is obvious that T < {{N^ - N)/2 + l). 

Note, that while estimating <3_d in this particular case we are taking into 
account the symmetry of the initial matrices. 

Thus the N-body system, according to the generalized S-tree method, 
is split into p-bounded subgroups, where p = {p^ ; p,i); (/ = 1, . . . , T) in the 
following way. We obtained matrix Du from matrix D with the help of 
some fixed p^ , afterwards we get matrix Du out of Dy^ using the described 
method. The last step of the algorithm is the distribution of the system into 
subgroups it via matrix for the fixed pi using the S-tree method. 

To get complete information of the distribution it is necessary to apply 
the algorithm for all the possible pairs {p^ ; p;), where p^ and pi accordingly 
are the components of the vectors pD and p. 

The time L, which is needed for the realization of this algorithm, is 
L < ((iV^ - N)/2 + 1)*+^ and for t << iV no difficulties can arise [|12|. 
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3 Parametrical approach 



As mentioned previously, the S-tree method is splitting the system into 
subgroups according to the boundness parameter p. 

The application of a parametrical approach to the S-tree method has 
two main goals 

1. To have a convenient way of representing output information of the 
S-tree algorithm with respect to individual properties of the particles. 

2. To have more complete and precise information about the substructure 
in the output information. 

Two possible versions of parametrical approach are considered. 

1. After obtaining the distribution of the system by subgroups (by the 
S-tree method), we introduce a parameter describing certain properties of 
particles independent of the distribution. 

Consider the set of N discrete points 

X = {xi, . . . ,XAr} 

and the function P, where 

P : X X X ^ R+. 

Let us introduce the following function 

B -.X ^ X ycR; 

B{xi) = {xi,bi). Wc denote {xi,bi) = x\^ for convenience, where bi is a 
parameter of P depending on the system under consideration. 

2. The suggested method is applicable if these constraints are satisfied 
simultaneously: 

1. k ^ N, where k is the number of subgroups; 

2. k ^ 1 and 3Zj, Ij, such that t;. > 1 and ti- > 1, where ti^ is the number 
of particles in the Ic-th subgroup, lc = l,...,k. 

For any fixed Xi the following sum is considered: 

where d is the number of Uj k subgroups and the summation is performed 

for all q = j, . . . ,k (accordingly, xj, . . . ,Xk forms the subgroup Uj^,,,^k and 

Xi 



6 



The sum 

Pi = ^P{Xi, Xq) 

(the summation is performed for all q = j, . . . , k, where Xi,Xj, . . . ,Xk form 
the subgroup C/j^...^fc) is considered as weh. The concept of parametric bound- 
ness is introduced in the Definition 3.1. 

Definition 3.1. If is more than any other Pf, than Xj is called a 
"parametrically real" (p.r) member of the subgroup; otherwise it is called 
"parametrically imaginary" (p.i). 

The final result of the output information in the parametric approxima- 
tion is presented in the following way. If .x,; is a p.r. member of subgroup 
Ul^, we denote Xi as x'*; if Xi is a p.i. member of subgroup Ui^, we denote Xi 
as , where Ij is the number of the subgroup Ui. with elements Xj, . . . ,Xk 
on which the function Pf{d = Ij) obtains its maximal value. 

4 Generalized parametric scheme 

Now we consider the corresponding generalizations for both versions of the 
parametrical approach. 

I. Consider the set of N points: X = {xi, . . . ,xn} and the function P, 
where 

P : X X X ^ R+. 

For the obtained distribution {via the S-tree method) we introduce a 
multidimensional parameter, each component of which shows a definite pe- 
culiarity of the element, independent of the distribution. 

Consider the following function 

B : X ^ X X R'',k= 1,2,... 

B{xi) = ixi,bi), 

where 6j = (bj, . . . ,b^). 

We introduce the notation {xi,bi) = x^' for convenience, where = 
1, . . . , is not a variable of P. 

In this particular case 6j is a multidimensional parameter, each com- 
ponent of which is normalized in accord with the physical content of the 
corresponding component and the system under investigation. 

II. The initial limitations of the proposed method fully correspond to the 
limitations of the parametrical approach. 
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We fix some value Xi and consider all possible sums Pf = T,P{xi,Xq), 
where d is the number of the subgroup Uj^^^^^k and the summation is taken 
for all the q = j, . . . , k (accordingly Xj, . . . ,Xk forms the group Uj^^^^^k and 

Xi 

It is obvious that the vector Pf = {Pf^, ■ ■ ■ ^~'')^ where k — z is the 
number of those subgroups for which the above mentioned limitations have 
been carried out. ^ 

We construct the vector P-^ from vector Pf, where Pf = {P- , . . . ,P^ ^^~^), 
the components of which are the components of the vector Pf situated in 
decreasing order. 

Let us consider the sum 

Pi = ^P{Xi, Xq) 

, (the summation is taken by all q = j, . . . ,k), where Xi,Xj, . . . , x^ forms the 
group C/i,„,,fc- 

The concept of the generalized parametric boundness is introduced in 
Definition 4.1. 

Definition 4.1 If Pf is greater then all components of the vector Pf, we 
say that Xi is "parametrically real" (p.r.) member of the group Ui^,_,^k, oth- 
erwise Xi is "parametrically imaginary" (p.i.) member of the group J7i,,,.,fc. 

The final representation of the output information in terms of generalized 
parametric approach is made in the following way: 

if Xi is a p.r. member of the group C//. , denote 

if Xi is a p.i. member of the group f7/. , denote / , where = 

. . .). Ij^ are numbers of those subgroups, where Pf^^ > Pf; r = 
1, . . . , k — z; dj^ = Ij^. 

5 Conclusion 

In this paper a parametric approach to the S-tree diagram method has been 
suggested. The proposed approach is useful for revealing the optimal distri- 
bution of the system under investigation |Q. 
The following developments were presented. 

• A generalization of the discrete S-tree method. 

• A generalization of the parametric approach to the S-tree method , in 
the case of multiparametric functions. 
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The methods described are more informative since they use all the infor- 
mation of the system in the input, while they enable the output to have a 
hierarchical structure depending not only on the kinematic and positional 
data, but also on the parameters having no direct influence on the dynamics 
of the system. 

For example, in the case of clusters of galaxies the algorithm will use data 
from the catalogues not only on the coordinates and redshifts of galaxies 
but also the morphological, color, spectral index, and other information 
self-consistently in the output information on the hierarchic structure of the 
systems. 

The hierarchical structure of N-body nonlinearly interacting particles is 
only one indication of their complex dynamics |13|,|14]. 
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